Composite Majorana Fermion Wavefunctions in Nanowires 
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We consider Majorana fermions (MFs) in quasi-one-dimensional nanowire systems containing nor- 
mal and superconducting sections where the topological phase based on Rashba spin orbit interaction 
can be tuned by magnetic fields. We derive explicit analytic solutions of the MF wavefunction in the 
weak and strong spin orbit interaction regimes. We find that the wavefunction for one single MF 
is a composite object formed by superpositions of different MF wavefunctions which have nearly 
disjoint supports in momentum space. These contributions are coming from the extrema of the 
spectrum, one centered around zero momentum and the other around the two Fermi points. As a 
result, the various MF wavefunctions have different localization lengths in real space and interfer- 
ence among them leads to pronounced oscillations of the MF probability density. For a transparent 
normal-superconducting junction we find that in the topological phase the MF leaks out from the 
superconducting into the normal section of the wire and is delocalized over the entire normal section, 
in agreement with recent numerical results by Chevallier et al.. 

PACS numbers: 73.63.Nm; 74.45.+C 
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I. INTRODUCTION 

Majorana fermions^ (MFs), being their own antipar- 
ticles, have attracted much attention in recent years in 
condensed matter physics^ - — . Besides being of funda- 
mental interest, these exotic quantum particles have the 
potential for being used in topological quantum comput- 
ing due to their non-Abelian statistics^ - —. There are a 
number of systems where to expect MFs, e.g. fractional 
quantum Hall system o 27 ' 28 , topological insulators 3 -^, op- 
tical lattices^, p-wave superconductors^, and especially 
nanowires with strong Rashba spin orbit interaction 7 - - — - 
the system of interest in this work. There are now sev- 
eral claims for experimental evidence of MFs in topolog- 
ical insulator a 29 ' 30 and, in particular, in semiconducting 
nanowires of the type considered here 3 -^ - — . 

As is well-known 7 - - — an s-wave superconductor 
brought into contact with a semiconducting nanowire 
with Rashba spin orbit interaction (SOI) induces effec- 
tive p-wave superconductivity that gives rise to MFs, one 
at each end of such a wire. Most studies have analyzed 
the corresponding model Hamiltonian by direct numeri- 
cal diagonalization, which provides exact solutions of the 
Schrodingcr equation for essentially all parameter values 
irrespective of their relative sizes. Less attention, how- 
ever, has been given to analytical approaches which can 
provide additional insights into the nature of MFs. As 
usual, this comes with a price: closed analytic expressions 
are hard to come by and can be obtained only in special 
limits. But since these limits turn out to include realistic 
parameter regimes such an approach is not a mere aca- 
demic exercise but worthwhile also from a physical point 
of view. 

Motivated by this, we focus in the present work on 
the spinor-wavefunction for MFs, and derive analytical 
expressions for various limiting cases, loosely character- 
ized as weak and strong SOI regimes. We find that these 
solutions are superpositions of states that come, in gen- 



eral, from different extremal points of the energy dis- 
persion, one centered around zero-momentum and the 
others around the Fermi points. Despite having nearly 
disjoint support in momentum-space, all such contribu- 
tions must be taken into account, in general, in order to 
satisfy the boundary conditions imposed on the spinor- 
wavefunctions in real space. As a consequence of this 
composite structure of the MF wavefunctions, there will 
be more than one localization length that characterizes 
a single MF. We will see throughout this work that the 
Schrodingcr equation for the systems under considera- 
tion allows, in principle, degenerate MF wavefunctions. 
However, this degeneracy gets completely removed by the 
boundary conditions considered here, and, consequently, 
there exists only one single MF wavefunction at a given 
end of the nanowire. The superposition also gives rise 
to interference effects that leads to pronounced oscilla- 
tions of the MF probability density in real space. Quite 
interestingly, the relative strengths of the different local- 
ization lengths as well as of the oscillation periods can be 
tuned by magnetic fields. 

If only a section of the wire is covered with a super- 
conductor, a normal-superconducting (NS) junction is 
formed. For this case, we find that the MF becomes de- 
localized over the entire normal section, while still local- 
ized in the superconducting section, as noted by several 
groups befor o 19 ' 34- — , and most recently studied in de- 
tail in a numerical study by Chevallier et al£&. Here, we 
will find analytical solutions for this problem, valid in the 
weak and strong SOI regime. Depending on the length 
of the normal section, the support of the MF wavefunc- 
tion is, again, centered at zero momentum or the Fermi 
momenta. Also similarly as before, different localization 
lengths and oscillation periods of the MF in the normal 
section occur, again tunable by magnetic fields. This 
could then provide an experimental signature for MFs, 
e.g. in a tunneling density of states measurement, where 
a signal that comes from a zero-mode MF will show os- 



2 



A, 



ical potential, the SOI term, 



FIG. 1. Nanowire (blue slab) of length L ac in the supercon- 
ducting regime with gap A sc induced via proximity effect by a 
bulk s-wave superconductor (not shown). A magnetic field B 
is applied along the nanowire in x-direction and perpendicu- 
larly to the Rashba SOI vector olr that points in z-direction. 



cillations along the normal section. 

The paper is organized as follows. In Sec. [II] we intro- 
duce the continuum model of a nanowire including SOI, 
magnetic field and induced superconductivity. The com- 
posite structure of MF in proximity-induced supercon- 
ducting wire is discussed in Sec. Mil for strong and weak 
SOI. In Sec. IIVI we investigate an NS junction and show 
how the type of MF wavefunction oscillates in space and 
depends on magnetic field. The final Sec. [V] contains our 
conclusions. Some technical details are referred to two 
Appendices. 



II. MODEL 

Following earlier work£~— i 19 i 38 ; our starting point is a 
semiconducting nanowire with Rashba SOI (see Fig. [1]) 
characterized by a SOI vector a R that points perpendic- 
ularly to the nanowire axis and defines the spin quantiza- 
tion direction z. In addition, a magnetic field B is applied 
along the nanowire in x-direction. We imagine that the 
nanowire (or a section of it) is in tunnel-contact with a 
conventional bulk s-wave superconductor which leads to 
proximity-induced superconductivity in the nanowire it- 
self, characterized by the induced s-wave gap A sc (see 
Fig. [T]). We refer to this part of the nanowire as to the 
superconducting section (or as to the nanowire being in 
the superconducting regime), in contrast to the 'normal' 
section of the nanowire that is not in contact with the 
superconductor and thus in the normal regime. 

We describe this nanowire system by a continuum 
model and our goal is to find the explicit wavefunctions 
for the MFs in the entire nanowire, including normal and 
superconducting section. For this, we need to introduce 
some basic definitions and briefly recall well-known re- 
sults about the spectrum. 

The Hamiltonian H° = H kin + H S01 + H z for the 
normal regime^ consists of the kinetic energy term 



H kin ^jdx *t(x) 



2m 
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H 



SOI 



-ia R J2 J dx yl(x){0 3 ) aa ,d x V a .{x), (2) 



where, again, the z-axis is chosen along a Rl and the Zee- 
man term corresponding to the magnetic field B along 
the nanowire (x-axis), 



= A; 



dx^iixXtn^^ix). (3) 



Here, ^\{x) is the creation operator of an electron at 
position x with spin a/2 = ±1/2 (along z-axis), and the 
Pauli matrices cti,2,3 act on the spin of the electron. The 
Zeeman energy is given by Az — gfisB/2, where g is the 
g- factor and \ib the Bohr magneton. It is convenient to 
introduce the corresponding Hamiltonian density H° , 

= *<(*)*>(,), 
H° = -h 2 dl/2m~ ii-ia R <7 3 d x + A z a u (4) 

which acts on the vector ip = (^^, ^i). The bulk spec- 
trum of H° (see Fig. [2£i) consists of two branches and is 
given by 



E° ± (k) 



2 1.2 



K 2 k 

2m 



f i±J(a R k) 2 + Ai 



(5) 



where k is a momentum along the nanowire. By open- 
ing a Zeeman gap 2A^, the magnetic field lifts the spin 
degeneracy at k = 0. The chemical potential [i is tuned 
inside this gap and set to zero. In this case, the Fermi 
wavevector is determined from E°_{kp) = and given by 



2kl 



4fcf D + k z , 



(6) 



where k so = ma R /H and kz = \/2Azm/h. 





where m is the (effective) electron mass and /i the chem- 



FIG. 2. Bulk spectrum for extended electron (solid lines) and 
hole (dashed lines) states in the normal (a) and in the super- 
conducting regime (b). (a) In the normal regime, a Zeeman 
gap 2Az is opened at k = 0, but the full spectrum is still 
gapless due to the propagating modes at the Fermi points 
±feir. (b) The proximity-induced superconductivity leads to 
the opening of a gap A e at the Fermi points ±Uf and modifies 
the topological gap A_ = A sc — Az at k — 0. 
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The nanowirc in the superconducting regime is de- 
scribed by the Hamiltonian H° + H sc , where the s-wave 
BCS Hamiltonian H sc couples states with opposite mo- 
menta and spinal, 

HSC = \Y, J dx A.c^.rCw-aW**' + h.c). (7) 

a. a' 

The proximity-induced superconductivity gap A sc is cho- 
sen to be real (thereby assuming that we can neglect the 
flux induced by the B-field, which is the case e.g. for 
InSb nanowires^). The spectrum of H° + H sc (see Fig. 
[2b) is then found to be 

El(k) = (^y + (a R kf + A% + Al (8) 

±2^A|AL+(^) 2 (A| + ( a ^). 

The 'topological' gap at k = is given by A_ = 
A sc — Az, and the closing of this gap marks the tran- 
sition between nontopological (A_ > 0) and topologi- 
cal (A_ < 0) phases^i^. In contrast, the gap at fcp, 
A e = 2\E-(Icf)\, is always nonzero (see Fig. [2p). 



III. MAJORANA FERMIONS IN THE 
SUPERCONDUCTING SECTION 

In this section we consider first the simpler case 
where the superconducting section extends over the en- 
tire nanowire from x = to x = L sc , see Fig. [TJ In 
the topological phase there is one MF bound state at 
each end of the nanowire^. In the physically interesting 
regime, these two MFs should be independent and have 
negligible spatial overlap. This justifies the consideration 
of a semi-infinite nanowire. In this work we focus on the 
MF at the left end, x = 0. 

We will consider two limiting regimes, namely strong 
(kp ~ 2k so ) and weak {hp — kz) SOI. In both regimes, 
the Hamiltonian can be linearized near the Fermi points 
and solved analytically. We show that the MF wavefunc- 
tion has support in fc-space from the exterior [k ~ ±fcp) 
and the interior (k ~ 0) branches of the spectrum, see 
Fig. [5Ji. If the system is in some intermediate regime 
of moderate SOI, the support of the MF wavefunction 
extends over all momenta from — kp to kp, and this case 
cannot be treated analytically in the linearization ap- 
proximation considered here. 

A. Regime of strong SOI and rotating frame 

The regime of strong SOI is defined by the condition 
that the SOI energy at the Fermi level is larger than the 
Zeeman splitting, Az *C maj^/H 2 (or kp w 2k so ), and 
larger than the proximity gap, A sc <C ma 2 R /h 2 . This 
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FIG. 3. (a) Band structure of a nanowire with strong SOI 
and in a uniform magnetic field B in the lab frame (see also 
Fig. [2]) for the normal section. States around k = belong to 
the interior branches and states around k — kF = 2k ao belong 
to the exterior branches, (b) The same band structure in the 
rotating frame. The rotating magnetic field B(x) given by 
Eq. H12|) couples Rf and Lj. that leads to the opening of the 
Zeeman gap 2Az but does not affect Lf and R±. 

allows us to treat the magnetic field and the proximity- 
induced superconductivity as small perturbations. 

The spectrum obtained in Eq. §5§ consists of two 
parabolas shifted by the SOI momentum k so = man/H 2 
and with a Zeeman gap opened at k = (see Fig. [5^,). 
In the strong SOI regime it is more convenient to work 
in the rotating frame, see Fig. \Sjp. For this we follow 
Ref. [H and make use of the following spin-dependent 
gauge transformation 

* a (x) = e-* k »**„(x), (9) 

where tilde refers to the rotating frame. The H S01 term 
is effectively eliminated (H S01 = 0) and the spectrum 
corresponding to H km consists of two parabolas centered 
at k = 0, one for spin up and one for spin down. Around 
the Fermi points, ±fc so , the spectrum can be linearized 
and the electron operators ^ a are expressed in terms of 
slowly-varying right (R a ) - and left (L a ) - movers, 

= Mx)e ik *° x + L a {x)e- ih >° x . (10) 
The kinetic energy term in the linearized model is 

H kin = ~ihv F J dx [Rl(x)d x R a (x) - Lt(x)d x L a (x)} 

(11) 

with Fermi velocity vp = a^/h. Here, we droped all fast 
oscillating terms, which is justified as long as £ 3> 2-k /k SOl 
where £ is a localization length of R a and L a (see below). 

In the rotating frame the B-field becomes helical, ro- 
tating in the plane perpendicular to the SOI vector an, 

B(x) = B[x cos(2fc so x) — y sin(2fc so x)]. (12) 

Here, x and y are unit vectors in x and y directions, re- 
spectively (see Fig. [T]). This leads to the Zeeman Hamil- 
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tonian of the form 



H z = A Z J dx (aOe 31 **"**-^*), 

~ A z I dx \F^{x)L^x) + L\(x)R^x)}, (13) 



where in the second line we used the linearization approx- 
imation and, again, droped all fast oscillating terms. We 
note that only R^(x) and L^(x) are coupled, which leads 
to opening of a gap, as shown in Fig. [3Jd. This is similar 
to the spin-selective Peierls mechanism discovered in Rcf. 
l39l where interaction effects strongly renormalize this gap 
(here, however, we shall ignore interaction effects). 

The superconductivity term [see Eq. (J7J] in the lin- 
earized model becomes 

H sc =i J dx A sc {R\(x)L±{x) - L^R^x) 

+ L^{x)R±{x) - R±{x)L t {x) + h.c). (14) 

We construct two vectors, </>^ = (Rf, L±, Rl, i|) and 
= (L^, R±, L|, which correspond to the exterior 
(k ~ 2k so ) and interior (k ~ 0) branches of the spectrum 
in the lab frame (see Fig. [3£i). The linearized Hamilto- 
nian H kin + H z + H sc reduces then to (l = i,e) 



= I j dx (^ l) (x)yH il) ^ l \x), (15) 



where the interior branches are described by 

U {l) = ~ihv F a 3 d x + A 2 criry 3 + A sc <r 2 i]2, (16) 

and the exterior ones by 

«W = ihv F a 3 d x + A sc a 2 r] 2 . (17) 

Here, the Pauli matrices 7/1,2,3 act on the electron-hole 
subspace. 

The energy spectrum is determined by the Schrodinger 
equation, 'H^tp E J(x) = E(pg(x), with boundary con- 
ditions to be imposed on the cigenfunctions (x) as 
discussed below. We introduce then the operator 7^ = 



/ dx tp® (x) ■ <f)^> (x) and see that it diagonalizes Eq. (fT5l) , 
i.e., = Ee e (JeV 

7^. Focusing now on the zero 



modes, we consider in particular 7''' 
it in a more convenient basis, 



'YeLo but express 



;(0 



(18) 



where \t = ^4,, ^j, and where fast oscillating 
terms were dropped. In this new basis ^P, we have rein- 
stalled the phase factors e ±tkpX (associated with R a and 

L a ) explicitly in the wavefunctions ^=o( x )' so * na * they 
are taken automatically into account when we impose the 
boundary conditions on ^^ =0 {x). 

The zero-energy operator 7") represents a MF, i.e., 
(7^))t = jy(0 5 if a nd only if the corresponding wavefunc- 

tion $e_ q {x) has the following form 



€Lo(*) 



(19) 



where the functions /, g are arbitrary up to normalization 
Jdx\S^ (x)\ 2 /2 = Jdx(\f(x)\ 2 + \g(x)\ 2 ) = 1, which, 
however, will be suppressed in the following. 

In infinite space (no boundary conditions), the spec- 
trum of the interior branches [see Eq. (|16l) ] is given 

by E<£> 



± 



Aj_, while the one for the 



(e) 



exterior branches [see Eq. (|17p] is given by E\ 
±y/ '(hvpn) 2 + Af , where k is taken from the Fermi point. 
Here, A± = A sc ± Az and A e = A sc . If Az and A sc 
become equal, the topological interior gap A_ is closed. 
In contrast, the exterior gap A e is not affected by the 
magnetic field. 

The only normalizable eigenstates of at zero en- 
ergy and at x > arc two evanescent modes com- 
ing from the interior branches, characterized by a de- 
cay wavevector h± = \A±\/an, and two evanescent 
modes coming from the exterior branches, characterized 
by a decay wavevector — A sc /a^>. The correspond- 
ing zero-energy eigenfunctions <Pe =0 (x) give the four ba- 
sis wavefunctions, exponentially decaying in the semi- 
infinite space x > 0, 
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f-i sgn(A_)e- ifc - :E \ 
i sgn(A_)e iks ° x 

„ — ik so x 



-k^'x 



I ^ e ik so x 



ik so x 



-i e 

Ak so x 



-i e 

Jtk 



ik so x 



I e ik so x ^ 

-i ( 7 



e + 



(20) 



(21) 



Here we should note that these four degenerate MF 
wavefunctions are not yet solutions of our problem: they 
do satisfy the Schrodinger equation but not yet the 
boundary conditions. Thus, we search now for a lin- 
ear combination of them, <Pm, such that the boundary 
conditions are satisfied. At the left end of the nanowire, 
the condition on the wavefunction is @m( x = 0) = 0. 
We assume here that the length of the nanowire L sc 
provides the largest scale, so we can neglect any inter- 
play between the two ends of the nanowire and treat 
them independently (see also below). The set of vectors 



W $W £(e) 



2 } is seen to be linearly independent 



in the nontopological phase at x = 0, thus it is impossible 
to satisfy the boundary conditions and no solution exists 
at zero energy. In contrast, in the topological phase, 
Az > A sc , the two vectors <P^ and are 'collinear' 
such that the boundary condition can be satisfied and the 
zero energy state is a MF given by <1>m = ^_ — ^f 1 in 
the rotating frame. Using Eq. ©, the MF wavefunction 
in the lab frame is then given by 









i e ikpX \ 








g—ikpx 


(1 






_^ e -ik F x 








e ik F x ^ 



(22) 



with fcp = 2k so . 

There are a few remarks in order. First, we see that the 
initial fourfold degeneracy of the MF has been completely 
removed by the boundary condition and we end up with 
one single non-degenerate MF wavefunction at the left 
end of the nanowire, x = (analogously for the right 
end, x = L sc ). This is reminiscent of a well-known fact in 
elementary quantum mechanics, where for spinlcss parti- 
cles in a one-dimensional box the degeneracy also gets re- 
moved by vanishing boundary conditions (whereas there 
is degeneracy for periodic boundary conditions) 42 This 
non-degeneracy of the MFs is a generic feature which will 
occur in all cases considered in this work, even in the 
presence of additional symmetries such as pseudo-time 
reversal invariance (see Sec. MI Bl below and footnote 

Second, we see that the MF wavefunction <I>m is a 
'composite' object that is a superposition of two MF 



I 

wavefunctions with (essentially) disjoint supports in k- 
space, one coming from the exterior ('P^) and one from 

interior ) branches of the spectrum, respectively. 
Note that the corresponding wavevectors are extrema 
of the particle-hole spectrum shown in Fig. ^jp. As a 
consequence, these two MF wavefunctions have differ- 
ent localization lengths in real space, 
£( e ' = l/k( e \ which are inverse proportional to the cor- 
responding gaps, |A_| and A sc . Which one of them de- 
termines the localization length of <P]\j depends on the 
ratio between Az and A sc . 

In particular, as the magnetic field is being increased 
from zero to the critical value B c = 2A sc /g/j,B, MFs 
emerge at each end of the nanowire^^. However, if the 
localization lengths £w of these emerging MFs arc com- 
parable to the length L sc of the nanowire, then these 
two MFs are hybridized into a subgap fermion of fi- 
nite energy (see App. lA"j) . This then implies that the 
MFs in a finite wire can only appear at a magnetic field 
B* ~ B c (l + Aan/ A SC L SC ) that is larger than the B c ob- 
tained for a semi-infinite nanowire. If the magnetic field 
is increased further, the main contribution to the MF 
bound state comes from the exterior branches. 

The composite structure of the MF wavefunction man- 
ifests itself in the probability density |<?m(^)| 2 along the 
nanowire. The density of a MF coming only from one 
of the branches, for example, <P_ , is just decaying ex- 
ponentially. In contrast, the density of the composite 
MF exhibits oscillations (see Fig. @|. These oscillations 
are due to interference and are most pronounced when 
the contributions of and <P_ to #m are similar, i.e., 
when both decay lengths, and S,^ , are close to each 
other. 

The approach of the rotating magnetic field allows us 
to understand the structure of composite MF wavefunc- 
tion. However, this approach is valid under the assump- 
tion that the SOI at the Fermi level is the largest energy 
scale. In order to explore the weak SOI regime, we come 
back to the full quadratic Hamiltonian in the next sub- 
section. 
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k so x 



h so x 



FIG. 4. The MF probability density \$ M (x)\ 2 , see Eq. {22}, 
for a nanowire in the strong SOI regime with h 2 A ac /ma 2 R — 
0.1. The decaying MF wavefunction $m oscillates with a 
period n/k ao . (a) In the topological phase but still near the 
transition, Az = 2A SC , the MF wavefunction undergoes many 
oscillations due to the interference between &_ and <P ^ . (b) 
Deep inside the topological phase, Az = 7A SC , the MF wave- 
function from the interior branches decays much faster than 
the one from the exterior branches. This leads to only a few 
oscillations of the density and a uniform decay with a decay 
wavevector fc' 6 ' away from the end of the nanowire. 



Weak SOI regime: Near the topological phase 
transition 



The regime of weak SOI is defined by the condition 
that the Zeeman splitting is much larger than the SOI 
energy at the Fermi level, Az 3> ma 2 R /h 2 [or kp as kz, 
see Eq. ©]. This allows us to treat the SOI as a pertur- 
bation. 

Around the Fermi points, ±/cp , the eigenstates of H° 
are found from the Schrodinger equation H°(—id x — > 
±k F )ip R/L = [see Eq. |@J] and given by 



R/L 



1 

71 



-i± 



1± «aa 



(23) 



where f> R ^ L denotes the eigenstates at k = ±kp. In 
Eq. (|23[) we kept only terms up to first order in k so /kp. 
As expected, (p R and (/?q are nearly 'aligned' along the 
magnetic field since Az ^S> a^kp. In the absence of SOI, 
lPq and ifiQ arc perfectly aligned along the x-axis and have 
the same spin, so they cannot be coupled by an ordinary 
s-wave superconductor. The SOI slightly tilts the spins in 
the orthogonal direction, which then allows the coupling 
between these states if the nanowire is brought into the 
proximity of an s-wave superconductor. 

The exterior branches can be treated in the linearized 
approximation similar to Sec. IIII Al 



X (x) = R{x)e lkFX + L{x)e 



—ikpx 



(24) 



where, again, R (L) annihilates a right (left) moving elec- 
tron. These operators are connected to spin-up (^-j-) and 
spin-down (^j.) electron operators as R(x) = ip R (x) ■ ip 
and L{x) = f L {x) ■ i/j, where ip = (^t'^-l) (with cor- 
responding support for right and left movers). Here, 
tp R / L (x) is given by Eq. (|23|) but where we allow now 
also for a slowly varying x-dependencc. 



In this approximation, we find 
H° = -ihv F [ dx [Rt(x)d x R(x) - L\x)d x L{x)], 



(25) 



where the Fermi velocity is given by vp ss h^/IAzjm. 
The proximity-induced superconductivity [see Eq. ([7])] is 
described in the linearized model as 

H sc = i J dx A sc (R(x)L(x) - L(x)R(x) + /i.e.), (26) 

where the strength of the proximity-induced effective p- 
wave superconductivity, A sc , is found from Eqs. ([7]). (|23p . 
and J2H), 



-r— = {<P0) ■ MT2<Po 



7. 

o "'so 

kp 



<C 1. 



(27) 



The suppression of A sc compared to A sc can be under- 
stood from the fact that two states with opposite mo- 
menta at the Fermi level have mostly parallel spins due 
to the strong magnetic field and they slightly deviate in 
the orthogonal direction due to the weak SOI, which then 
leads to a suppression of A sc by a factor k so /kp. 

Again, introducing a vector <fi^{x) = [R,L,R) ,L^), 
we represent the linearized Hamiltonian = H° + H sc 
as 



#(<o = I J dx(^n^U {e \ 

= ~ihv F T 3 d x + A sc T 2 r?2, 



(28) 



where the Pauli matrices Ti,2,3 act on the right/left- mover 
subspace. 

The spectrum around the Fermi points in infinite space 
(no boundary conditions) follows from the Schrodinger 

equation, %^ip^ = E^(p^ e \ and is given by = 

± y (Hvfk) 2 + Ag, where the momentum k is again taken 
from the Fermi points. Here, 2A e = 2A SC is the gap 
induced by superconductivity. 

The zero-energy solutions that are normalizable for 
x > are two evanescent modes with wavevector 
fc( e ) = A sc /hvp determining the localization length. 

These solutions can be written explicitly as (p± = 

(1, — i, l,i)e~* " )x and (p^ = (— i, 1, i, l)e~^ Ce x . Repeat- 
ing the procedure that led us to Eq. (fl~8|) in Sec. IIII Al 
we can introduce a new MF operator, 



7 = / dx <P E =o(x) ■ <J(x). 



(29) 



Here \P = (^f , ^4., ^j). The two corresponding wave- 
functions arc written as 



<I> 



(e) 



/;(*) 

\-if,i*V 



-fe (e) x 



(30) 
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where 

fi(x) = i(l + k so /k F )e lkFX - (1 - k so /k F )e- ikFX , 
f 2 {x) = i(l - k ao /k F )e- ikpX - (1 + k so /k F y kFX . (31) 



The effect of the SOI on the states around fc = 
is negligible near the topological phase transition if 
h 2 \ A_ |/2ma|j -C 1. Therefore, the eigenstates for weak 
and strong SOIs are the same in first order in SOI. This 



means that we are allowed to take S_ and given by 
Eq. (|20p and transform them back into the lab frame, 




After we found four basis wavefunctions 
,$[ e \^ 2 e) }, we should impose the bound- 
ary conditions on their linear combination <P(x). The 
wavefunction <P(x) should vanish at the boundary x = 0. 
One can see that if we neglect the corrections to the 
wavefunctions <P^ and <P^ coming from SOI, then we 
are able to satisfy the boundary conditions. This is a 
consequence of the fact that in the absence of SOI both 
states at the Fermi level have the same spin, so that 
the functions <p\ e \ i = 1,2, effectively become spinless 
objects and the MF always exists and arises only from 
the exterior branches. For the complete treatment, 
however, we should also consider contributions from the 
interior branches. This will be addressed next. 



The set of wavefunctions {<£ 



W *« rf,^) 



comes linearly dependent in the topological regime A_ < 
and the MF wavefunction is given by 

$ M = ( 1 - M 4 e) - ( 1 - M 4 e) -4^#« . (33) 
V k F I \ k F I k F 



As in the regime of strong SOI (see Subsec. MI Al) . 
the MF wavefunction has its support around wavevec- 
tors k = (interior branches) and k = ±k F (exterior 
branches). However, in contrast to the previous case, the 
contribution of the interior branches is suppressed by the 
small parameter k so /k F , thus the exterior branches con- 
tribute most to the MF wavefunction. At the same time 
we note that the localization length of MFs is determined 
by the smallest gap in the system. Near the topological 
phase transition^, which corresponds to the closing of 
the topological (interior) gap, the interior branches deter- 
mine the localization length as long as k^ > k_ . If the 
magnetic field is increased further, the gap in the system 
is given by the exterior gap, 2A e = 2A SC oc 1/VB [sec 
Eq. (|27p]. The localization length is increasing as oc \J~B. 
As soon as it is comparable to the nanowire length L ac , 
the wavefunction of the two MFs at opposite ends over- 
lap, and the two zero-energy MF levels are split into one 
subgap fcrmion of finite energy. 



\ I^m(*)| 2 




^V-^ k F i 
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FIG. 5. The MF probability density \$m(x)\ 2 for a nanowire 
in the weak SOI regime ( ma^/h 2 A sc = 0.2 and Az = 5A SC ) 
oscillates with period tt /k F due to interference between right- 
and left-moving contributions [see Eq. (|30[) ]. The decay 
length is given by l/fc (e) . 



In the weak SOI regime and sufficiently far away from 
the topological transition point, Az > A sc (l + k so /k F ), 
so that the gap is determined by the exterior branches 
only, we can work in the simplified models given by 

The explicit MF wavefunction can be 
dS3J), 



[see Eq. ([28])]. 
found from Eq. 



jgiTT/4 

e j7r/4 



sin(fci?a;)e 



(34) 



Again, we note that this wavefunction describes a MF 
with the spin of both, the electron and the hole, point- 
ing in x-direction, again, up to corrections of order 



§>\ '} be- of k so /k F . The MF probability density \$ M {x){' 



sin (k F x)e 



decays oscillating with a period half 



the Fermi wavelength, \ F /2 = ir/k F (see Fig. [5J. 

In passing we remark that given in Eq. (|28p be- 
longs to the topological class Dili according to the clas- 
sification scheme of Ref. and supports MFs in ID, in 
agreement with our result Eq. (|34t 41 . 



OtR 




A.. 



-L r 



N 



FIG. 6. NS junction of a nanowire. The right section (blue) 
of the nanowire from x = to x = L BC is brought into contact 
with a bulk s-wave superconductor (not shown) that induces a 
gap A ac in the nanowire via proximity effect. The left section 
(green) of the nanowire from x — —L n to x = is in the 
normal regime. A magnetic field B is applied along the entire 
nanowire in i-direction and perpendicularly to the Rashba 
SOI vector olr, which points in z-direction. 
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IV. MAJORANA FERMIONS IN NS 
JUNCTIONS 

In this section we consider a nanowirc containing a 
normal-superconducting (NS) junction where the right 
part is in the superconducting and the left part in the 
normal regime, see Fig. [5] The junction is assumed to 
be fully transparent. We will show that the MF wave- 
function leaks out of the superconducting section and 
leads to a new MF bound state that extends over the 
entire normal section. We note that this bound state is 
different from Andreev bound states 4 ^ known to occur in 
NS junction systems. Indeed, the existence of the latter 
at zero energy would be accidental in the presence of a 
magnetic field since they move away from the Fermi level 
if the magnetic field is varied. Further, the MFs found 
in this section always exist in the topological regime and 
are not sensitive to the length L n of the normal section, 
in stark contrast to Andreev bound states that move in 
energy as function of L „ 43 i 44 . 

We continue to work with the formalism developed in 
Sec. |TlT] [see Eqs. (fT8|) and (|29)) ] and represent 7 in the 
basis of electron/hole spin- up/spin-down operators ^(2) 
in terms of a four-component vector <P(x) on which we 
impose the boundary conditions. As before, the length 
of the superconducting part of the nanowirc L sc is as- 
sumed to be much larger than any decay length given 
by k_ , k+\ or k^ (fc( e '). This assumption allows us to 
treat the nanowire again as semi-infinite with no bound- 
ary conditions at x = L sc . In contrast, the normal sec- 



tion, x € [— L n ,0], is finite. Thus, at x = we invoke 
continuity of the wavefunctions and their derivatives and 
at x — —L n we impose vanishing boundary conditions, 

${x = 0") = $(x = 0+), (35) 
d x ${x = 0") = d x $(x = 0+), (36) 
$(x = -L n ) = 0. (37) 

The analytical form of the functions <P(x) can be found 
in two regimes, again in the weak and strong SOI limits, 
which we address now in turn. 

A. NS junction in the strong SOI regime 

As before, the most convenient way to treat the strong 
SOI regime is to work in the rotating frame. In Sec. 
IIII Al we already found the four basis wavefunctions at 
zero energy 

(${ sc) , $ { 2 SC) , $ { 3 SC) , $[ sc) ) = , , $[ e) , $ { 2 e) ) (38) 

in the superconducting section, for x > [see Eqs. ([20)) 
and ([21])]. 

The eigenfunctions for the normal section can be 
found from the_ linearized Hamiltonians for the inte- 
rior branches, [see Eq. (|16p], and for the exterior 
branches, HS e ^ [see Eq. (fT?)) ]. with A sc = 0. The exterior 
branches are not gapped leading to the four propagating 
modes (see Fig. ^jp) with wavefunctions given by 



(n) 



( —% e 

e~ 
i e 
e 



k so x 
k so x 
ik so x 
ik so x 



/ „ifc 

V 



e""'° x \ 

— ik so x 



( I 



i e 



ik so x 



(n) 



-i e %ks ' 

Ak so x 



\ 

c 

/ 



(n) 



—i e 



y i e 



ik so x 



(39) 



where we choose to represent the wavefunctions in form 
of MFs, guided by our expectation that the final solution 



is also a MF. The interior branches are gapped by the 
magnetic field (see Fig. [3Jd) and the four corresponding 
wavefunctions describing evanescent modes are given by 



(«) 



-ik so x 



^ n) x 



-ik s 



Z """ x \ 
„ik so x 



% e 

„ik so x 



(40) 



-ik BO x 



(n) 



-1 e s ° 
\ e - ik "" x ) 



( e - ik >° x \ 



(n) 



-1 e 



ik* 



\i e - ik >° x ) 



-k (n) {L n +x) 



(41) 



with fcW = Az/ctR. The modes decay from their maximum at x = to zero for x — > —00, while decay 
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we search for their linear combination, 



FIG. 7. The MF probability density |<£ A /(a;)| 2 in an NS junc- 
tion for a nanowire in the strong SOI regime (ft 2 A sc /ma:f{ = 
0.06). The normal section of length L n is long compared to 
the decay length, i.e. k^L n 2> 1. The MF wavefunction ex- 
tends over the entire normal section (green) and decays expo- 
nentially inside the superconducting section (white). The os- 
cillations with period n/k so result from interference between 
the three components r> °f in wea k magnetic fields 
(Az = 1.5A SC , blue full line), fc' n 'L n ~ 1, and the oscillations 
extend over the entire normal section. In contrast, in strong 
magnetic fields (Az = 7A SC , red dashed line), the oscillations 
are strongly suppressed. Inset: |<?m(:e)| 2 as function of x for 
a short normal section, k^ n 'L n <C 1. Similarly to the case 
of the superconducting wire (see Fig. Q, \<Pm(x)\ 2 decays 
oscillating. 



&m{x) 



V 4 a $ {sc) 



V 8 b <P {n) 



x > 

-L n < x < 0, 



(42) 



such that the boundary conditions (|35l) - ([37|) are satis- 
fied. This is, in general, possible only in the topological 
phase. However, we also find solutions in the nontopolog- 
ical phase, where these solutions exist only if some special 
relations between the parameters Az, cur, A sc , and L n 
are satisfied. This allows us to identify them as Andrccv 
bound states in an NS junction. Since they are not of 
interest here, we focus on the solutions in the topological 
phase only. 

It is worth of pointing out that due to the internal sym- 
metry of the MF wavefunctions, five of the coefficients are 
readily seen to vanish, namely 



a 2 = h = b 4 



0. 



(43) 



The exact analytical solution is given in App. [B] and 
used for the plot in Fig. [7] Here, we only discuss the two 
limiting cases of long and short normal sections L n . 
First, we consider L„ 3> l/k^ n \ allowing us to neglect 



the terms and at x 



-L n , and and 



from their maximum at x 



-L n to zero for 



+oo. 



at x — 0. In this case, the sum in Eq. (|4"2"j) is determined 
by the coefficients 

oi -4-0, b 6 -> 0, b 7 = -1, (44) 
a.3 = h = cos(2k so L n ), a 4 = b 2 = - s'm(2k so L„), 



After having introduced the basis consisting of 12 MF 
wavefunctions given by Eqs. (|2"0]). (j2Tj) , ([55]), and PT|). 



leading to the solution of the form 



<4 n) {2L n + x) e-^° )x . 
$i n \2L n + x)-$ { 3 n \- 



3 -feW(x+L„) 



x > 

-L n < x < 0. 



(45) 



Thus, we see that the MF wavefunction in the lab frame, 
&m(x), decays monotonically in the superconducting 
section while it oscillates in the normal one (see Fig. 
[?])• In weak magnetic fields the MF probability den- 
sity |<?m(^)| 2 oscillates over the entire normal section, 
in contrast to the near absence of oscillations in strong 
magnetic fields. 

We note that a long normal section serves as a 'mo- 
mentum filter'. As shown in Subscc. IIII A[ a MF has 
equal support from the exterior and interior branches [see 
Eq. (|22j) ] if the entire nanowire is in the superconducting 
regime. In contrast to that, if a significant portion of 
the nanowire is in the normal regime, the MF has sup- 
port mostly from the exterior branches with momenta 
k ~ ±kp. The contributions from the interior branches 
with momenta k ~ arc negligibly small, ai — > and 



a 2 = 0. This behavior can be understood in terms of 
momentum mismatch: the normal section does not have 
propagating modes with k — (in the lab frame). Thus, 
while such k = modes exist in the superconducting 
section, they cannot propagate into the normal section. 

Second, we consider the opposite limit L„ -C l/k^ n \ 

— (n) 

Here, we can treat the decaying solutions, <&)■ , j = 
5,6,7,8, as being constant over L n . The MF wavefunc- 
tion $m(x) is constructed from seven basis MF wave- 
functions with the same coefficients as in Eq. (|4"4")) with 
the only difference that now a\ = — 1. For short normal 
sections, the form of the MF probability density |<?m(^)| 2 
is very similar to the one of a superconducting nanowire 
(compare inset of Fig. [7] with Fig. 0^.). Again, the inter- 
ference between , , and leads to oscillations 
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in the superconducting section. 

In both limits of short and long normal sections, 
|^m(^)| 2 has it maximum in the normal section while 
it decays in the superconducting one. This opens the 
possibility of measuring the presence of a MF state spec- 
troscopically in the normal section. The amplitude and 
period of oscillations of the MF probability density is 
sensitive to magnetic fields and to the nanowire length, 
controlled e.g. by an infinite barrier on the left end. 
Moreover, by shifting such a barrier via gates, we can 



change the type of MF from $Y to (P^' [ see Eqs 



rAe) 



(|42j) . and (|44|) ] . This amounts to change a given MF state 
from a 'real part' type, tp + ijj', to an 'imaginary part' 
type, i(ip - V-' 1 )- 



B. NS junction in the weak SOI regime 

As before, we first identify basis wavefunctions in 
the superconducting section and in the normal section. 
Then, we search for a linear combination of them such 
that the boundary conditions given by Eqs. (|55|) - (|37p arc 



satisfied. 

As shown in Sec. IIIIB| the MF wavcf unction #m 
has predominantly support from the exterior branches. 
The correction to the MF wavefunction from the interior 
branches is suppressed by a factor k so /kF [see Eq. (|33|) ]. 
If we focus on the regime away from the topological phase 
transition where the exterior gap is smaller than the inte- 
rior one, then, as in Eq. (|50"|) . the MF wavefunction can 
be constructed to first order in k so /kp from the exterior 



wavefunctions <P 



3=1,2 



alone. 



fji 



(sc) 



( 9 3 {x) \ 
\-igj(x)J 

- in / 4 sm(k F x), 52 = e-^cos^Fx). (46) 



The propagating electron modes ip R and <p L of Eq. 



([23]) in the normal section described by H° (without H s 
were considered before [see Eq. (j2"9")) ] and given by 



(n) 



— ikp-x 
e -ik F x 





f e ik F x \ 




/ ie -ik F x 




/ ie ik F x \ 


, 4 l) - 


^ik F x 

g — ik F x 
\ e -ik F xj 


, *w = 


_ ie ^k F x 

—ie ifCFX 
1 ie ik F x j 


' ^4 


-ie lkFX 

_i e -ik F x 
1 ie -ik F x 



(47) 



The ansatz for the wavefunction in both sections is 



(sc) 



2^3=1 "^3 



X > 

-L n <x<0. 



(48) 



The coefficients a 3 and bj can then be found from the 



boundary conditions (f35 j) -([37 j) . 

a% = 2[k~( e > sm(kFL n ) + kF cos(fc_pL„)], 

a 2 = 2kF sin(fcF^n), 

b\ = —64 = kp cos(kpL n — 7r/4), 

62 = -b 3 = -k F cos(k F L n + tt/4), (49) 

leading finally to the MF wavefunction of the form 



«MaO = f(x) 



.g-iTT/4 



(50) 



where 



kp sm(kp[x + L n ]), 



-L n < x < 



-k e x 



[kp sin(kp[x + L r 



-k e sm(kpx) sin(fci?L n )], x > 0. 



(51) 



The corresponding MF probability density |$m(2:)| 2 is 
shown in Fig. |U The MF wavefunction extends over the 
entire normal section. In this section, we considered as 

(n) 

basis functions only propagating modes, <Pj , leading to 
a purely oscillatory solution with the period given by half 
the Fermi wavelength Xp/2 = n/kp. In contrast, in the 
superconducting section, the MF wavefunction decays on 
a short distance. In other words, the MF is mostly de- 




FIG. 8. The MF probability density \$ M (x)\ 2 in an NS junc- 
tion for a nanowire in the weak SOI regime (ma^/h 2 A sc = 
0.2 and Az = 5A SC ). The MF wavefunction extends over the 
entire normal section of length L n (green) and decays rapidly 
inside the superconducting section (white). 



11 



localized over the entire normal section and is strongly 
localized in the superconducting section, in agreement 
with recent numerical results^. This might simplify the 
detection of MFs by local density measurements since the 
normal section is freely accessible to tunnel contacts, in 
contrast to the superconducting section which needs to 
be covered by a bulk s-wave superconductor. 



V. CONCLUSIONS 

In this work, we have focused on the wavefunction 
properties of Majorana fermions occurring in supercon- 
ducting nanowircs and in nanowires with an NS junction. 
The superconducting phase is effectively p-wave and is 
based on an interplay of s-wave proximity effect, spin 
orbit interaction, and magnetic fields. We have derived 
explicit results for the MF wavefunctions in the regime 
of strong and weak SOI and shown that the wavefunc- 
tions are composite objects, being superpositions of con- 
tributions coming from the interior (around k — 0) and 
exterior (around ±fep) branches of the spectrum in mo- 
mentum space. While the underlying Hamiltonians con- 
sidered in this work allow degenerate MF wavefunctions, 
the boundary conditions at hand completely lift this de- 
generacy and we are left with only one single MF state 
at a given end of the nanowire (i.e. in total there are two 
MF states for the entire nanowire). 

In the strong SOI regime of a superconducting 
nanowire both branches contribute equally. However, the 
decay length of the MF is determined by the branch that 
also defines the smallest gap in the system. Moreover, the 
oscillations in the MF probability density with period of 
the Fermi wavelength are decaying on the scale given by 
the largest gap in the system. In the weak SOI regime, 
the exterior branches mostly contribute to the MF wave- 



function. The contributions of the interior branch is sup- 
pressed by the small factor k so /kp <C 1 and only close to 
the topological phase transition this branch determines 
the localization length of the MF. The interference be- 
tween modes from kp and —kp leads to oscillations of 
the probability density of the MF with an exponentially 
decaying envelope. 

For a nanowire with an NS junction we find that the 
MF wavefunction becomes delocalized over the entire 
normal section, while still being localized in the super- 
conducting section, in agreement with recent numerical 
results^ Again, we obtain analytical results for the weak 
and strong SOI regimes. Depending on the length of 
the normal section, the support of the MF wavefunction 
is centered at zero momentum or at the Fermi points. 
Again, we find different localization lengths and oscilla- 
tion periods of the MF in the normal section that are 
tunable by magnetic fields. Based on this insight, we ex- 
pect that in a tunneling density of states measurement 
the tunneling current at zero bias exhibits oscillations as 
a function of position along the normal section due to the 
presence of the MF in the normal section. 

Finally, we remark that in this work we have focused 
on single particle properties and ignored, in particular, 
interaction effects. It would be interesting to extend the 
present analysis to interacting Luttinger liquids, in par- 
ticular for the SOI nanowire with an NS junction, com- 
bining the approaches developed in Refs. Q, [H, HI, and 

m 
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Appendix A: Finite nanowire 

In this Appendix we address the problem of a finite 
superconducting section of length L sc . In this case, 

the decaying modes <P = , ^+\ ^i^, @i } Wltn 
maximum at x = are given by Eqs. (|20]l and pip . 
Now we should also take into account the four evanes- 
cent modes with maximum at x — L* r . These modes, 



<I> 



(e) 



$2^}i are f° un( i from Eqs. 



l"6|) . (fTTf and are similar by their structure to S , 



(») 



(e) 



-i sgn(A_)e ifcsoa: 

ik so x\ 



-i e 



-ik so x 



i e 



-ik so x 



k™(.x-L ac ) 



g — ik so x \ 

i e lks ° x 

■i e ~ lk soX j 

^ik so x \ 
ik so x 



-i e 
e 



ik so x 
k. 



-kf\x-L sc ) 



-k {c) (x~L sc ) 



(Al) 



We construct a 8 x 4 matrix uj(x) from the eight ba- 
sis wavefunctions. The zero-energy solution can then be 
compactly written as 

*(x) - a • $ (x) + b • $ Lsc {x) = u(x) Q , (A2) 

where a = (01,02,03,04) and b = (b±, 62, 63, 64) are co- 
efficients that should be determined from the boundary 
conditions, 

$(x = 0, L sc ) = 0, (A3) 



which can be rewritten as a matrix equation in terms of 
a 8 x 8 matrix L7, 

W = (2(0) 2(L. C )) (b)=°- ( A4 ) 

The determinant of the matrix £1 is nonzero, so the 
solution of the matrix equation is unique and trivial, 
(a, b) = 0. This means that, strictly speaking, MFs 
cannot emerge in a finite-size nanowire. MFs exist only 
under the assumption that the overlap of the two MF 
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wavcfunctions (localized at each end of the nanowire and 
derived in a semi-infinite nanowire model) can be ne- 
glected. Otherwise, the two MFs are hybridized into a 
subgap fcrmion of finite energy. 

Appendix B: Exact solution in strong SOI regime 

Here, we present the exact solution for the MF wave- 
function $(x) composed of seven different basis MFs 

I 



wavefunctions [see Eq. (|42|) ] and satisfmg the boundary 
conditions ([35 ]) - ([3T]t . We find 



a 2 = bi = 64 = h = h = 0, 

ai = Ak 2 a + Ak so k (e) cosh(fc ( " ) L) sin(2fc SOJ L n ) - 2fc (e) cos(2fc so L„)[fc (e) cosh(fc ( " ) L„) + fc (n) sinh(fc ( " ) i„)], 
a 3 = e - fc(n)L -[e ftWi "(feW[fcW _ 2 fc«] + jfeW[jfcC») - fc«]) - 4k 2 so e 2ki " ] L " cos(2k so L n ) 
- 2k so {e 2kin)L -k {e) + fc W - fc (Tl) )sin(2fc SOJ L„)], 

a 4 = -e- fcl " >L "[2e fc( ™ >L "fc so (fc(") - k ( l ] - k^) + 2k so (e 2k{n) L "k^ + fe<? - k {n) ) cos(2fc so L„) 

_ e 2fc<-)£ B4fc 2 o sin ( 2 fc so i„)] ; 

fo 2 = -e- fc( " )i "(-2e fc( " >i "A :so fc( e ) + 2k so {k ( l ] - *W) cos(2fc so L n ) - 4fc 2 e 2fe< " ,L " sin(2fc so L„)), 
63 = e - fc( " ,L "(e fc( " ,L "fc( e )[fc^ - - 4fc 2 e 2fe( " )L " cos(2fc so L„) + 2fc so (fc<") - fc W ) sin(2fc ao L n )), 
b 6 = -2k so (k^ - k ( l ] ) - e fc( " >i "fc^[2fc so cos(2fc so L„) + (k^ - fc W ) sin(2fc so i„)], 

fe 7 = 4fc 2 e fc( " )L ". (Bl) 



